Scalar field self-force effects on orbits about a Schwarzschild black hole 
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For a particle of mass fi and scalar charge q, we compute the effects of the scalar field self-force 
upon circular orbits, upon slightly eccentric orbits and upon the innermost stable circular orbit of 
a Schwarzschild black hole of mass m. For circular orbits the self force is outward and causes the 
angular frequency at a given radius to decrease. For slightly eccentric orbits the self force decreases 
the rate of the precession of the orbit. The effect of the self force moves the radius of the innermost 
stable circular orbit inward by 0.122701 x q 2 /fx, and it increases the angular frequency of the ISCO 
by the fraction 0.0291657 x q 2 /fim. 
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I. INTRODUCTION 

We consider a small mass /i with a scalar charge q 
which orbits a Schwarzschild black hole. The interac- 
tion of the particle with its own field yields the self force 
on the particle. We use a perturbation analysis to find 
the effects of the scalar-field self-force upon circular or- 
bits, upon slightly eccentric orbits and upon the loca- 
tion of the innermost stable circular orbit, the ISCO, of 
the Schwarzschild geometry. We view this project as a 
warm-up for the more interesting gravitational problem 
which must deal with more complicated field equations 
and gauge issues [1-3]. 

In general relativity, a particle of infinitesimal mass 
moves through spacetime along a geodesic. If the parti- 
cle has a small but finite mass fj, then its worldline devi- 
ates from a geodesic of the background spacetime by an 
amount proportional to fi. This deviation is said to result 
from the "self force" of the particle's own gravitational 
field acting upon itself. 

Newtonian gravity presents an elementary example of 
a self- force effect [4]. A small mass /i in a circular orbit 
of radius r about a more massive companion m has an 
angular frequency fl Q given by 

= 3n " , 2 - (1) 

When fi is infinitesimal, the large mass m does not move, 
the radius of the orbit r D is equal to the separation be- 
tween the masses and fl a 2 — m/r 3 . However when fx 
is finite but still small, both masses orbit their common 
center of mass with a separation of r Q (l + /u/ra), and the 
angular frequency is as given in Eq. (1). The finite fj, 
influences the motion of m, which influences the gravita- 
tional field within which /i moves. This back action of /i 
upon its own motion is the hallmark of a self force, and 
the n dependence of Eq. (1) is properly described as a 
Newtonian self force effect. 

A thorough understanding of gravitational waves de- 
tected by the Laser Interferometer Space Antenna [35] 
will require clear theoretical predictions of possible grav- 



itational wave forms which result from a small stellar- 
mass object orbiting a jumbo sized black hole; these wave 
forms must include self-force effects. 

In studying the gravitational self force, one considers 
the particle's gravitational field to be a small perturba- 
tion h a b of the background metric g a b- For an object of 
very small size, the motion ought to be independent of 
the particle's structure, and one is inclined to take the 
limit of a point particle. However, in that limit, h a b di- 
verges precisely at the particle, and the concept of the 
self force might appear to be ill-defined. 

Dirac [5] studied the electromagnetic version of this 
problem in flat spacetime and discovered that the part 
of the actual, retarded electromagnetic field which is sin- 
gular and yet exerts no force on the particle itself is, in 
a local approximation, the Coulomb field and could be 
identified as the average of the retarded and advanced 
electromagnetic fields. The remainder, half of the dif- 
ference between the retarded and advanced fields, is a 
vacuum solution of Maxwell's equations and accounts en- 
tirely for the particle's self force. 

The self force includes the radiative reaction force or 
radiation reaction [6]. For a particle with electrical charge 
q and acceleration a, the Abraham-Lorentz force de- 
scribes the response of a particle to its own radiation and 
is proportional to q 2 da/dt. The factor of q 2 results from 
the charge q interacting with its own electric field, which 
is also proportional to q. Similarly the gravitational radi- 
ation reaction force on a small mass fj, is proportional to 
li 2 . Other parts of the self force are not directly related 
to radiation but are properly described as the particle 
interacting with its own field and are also proportional 
to q 2 or fi 2 . 

A curved-spacetime generalization of Dirac's approach 
is now available [7-11]. For gravitation, an expansion 
about the position of the particle describes the singular 
"S" part of the field h^ b which exerts no force on the 
particle and is a local solution of the perturbed Einstein 
equations with the particle as the source. The remainder 
"R" part of the field h^ b = h a b — h^ b is, locally, a source- 
free solution of the perturbed Einstein equations, with 
the combined metric g a b + h^ b being a vacuum solution 
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of the Einstein equations through first order in /i R fc . The 
effect of the self force has the particle moving along a 
geodesic of the vacuum geometry g a b + h^ b . 

A caveat remains: in curved spacetime the "S" and 
"R" fields cannot be described in terms of the advanced 
and retarded fields. However, the mode-sum regulariza- 
tion procedure pioneered by Barack and Ori [12-21] is of 
use for background geometries amenable to the decompo- 
sition of fields in terms of scalar, vector and tensor har- 
monics. This procedure has been applied to self- force cal- 
culations for the Schwarzschild geometry involving both 
scalar [19, 22, 23] and gravitational [24] fields. In gen- 
eral terms, the "S" field is singular at the particle, but 
each of its spherical harmonic components is finite; these 
components are the "regularization parameters". The 
regular "R" field is determined by subtracting each £, m 
mode of the "S" field, from the corresponding £, m mode 
of the actual retarded field. The sum over modes of the 
difference between the retarded and "S" fields provides 
the "R" field which governs the self-force effects on the 
particle. 

In this manuscript we treat the self force from the 
scalar field $ in a pcrturbative manner. At zeroth or- 
der in 'J, a geodesic for the particle is chosen and this 
determines the particle's singular field \& s in the neigh- 
borhood of the geodesic. Then the actual scalar field 'S, 
with appropriate boundary conditions, is found every- 
where. The difference of the actual field and the singular 
"S" field provides the regular remainder ^ R = ^ — 
Finally we determine the first order in 5* self-force effects 
of the scalar field <3/ R acting back on the particle and 
changing the worldline away from the original geodesic. 
The small effects of ^ R appear as an acceleration of the 
worldline of order g 2 //im 2 . We determine the effects of 
the scalar field self-force upon the angular frequency, en- 
ergy and angular momentum of a general circular orbit, 
upon the rate of precession of a slightly eccentric orbit 
and upon the location and angular frequency of the in- 
nermost stable circular orbit. 

In Section II we give an overview of the self force and 
how it affects the location of the ISCO. 

In Section III we describe how a generic scalar field 
affects the worldline of a particle in the Schwarzschild 
geometry. 

In Section IV we consider scalar field self-force effects 
on circular orbits, and numerically calculate the corre- 
sponding changes of energy, angular momentum and an- 
gular frequency. In particular the self-force effect on the 
right hand side of Eq. (4), below, is derived. We also give 
the results of a numerical computation of , 3> R evaluated 
at the location of the particle, which changes the effective 
inertial mass \i of the particle. 

In Section V we consider the scalar-field self-force ef- 
fects on slightly eccentric orbits of the Schwarzschild ge- 
ometry. Specifically, we find the self-force correction to 
the rate of the precession of the orbit and to the right 
hand side of Eq. (5), below. In Section VI, this analy- 
sis is applied to an orbit near r — 6m to obtain the self 



force effect upon the radius and angular frequency of the 
ISCO. 

The discussion of Section VII distinguishes between 
the scalar and gravitational self force effects. We do not 
attempt to generalize our limited results to the case of 
the gravitational self-force. 

In Appendix A we give some details of the spheri- 
cal harmonic decomposition of a source moving along a 
slightly eccentric orbit. These details are required for 
the numerical calculation of the actual field for a slightly 
eccentric orbit. 

In Appendix B we describe the regularization param- 
eters which are required for the self-force analysis of 
slightly eccentric circular orbits of the Schwarzschild ge- 
ometry. These are used in the computation of \I/ R de- 
scribed in Sections IV- VI. 

We use Schwarzschild coordinates (t, r, 8, </>) on the 
spacetime manifold. The position of the particle in these 
coordinates is (T, R, 6,$), and s measures the proper 
time along a worldline. For a general worldline, the an- 
gular frequency is £1^ = d^/dT. The subscript is re- 
served exclusively for quantities related to circular orbits: 
for a circular orbit, the radius and angular frequency are 
r Q and f2 > respectively, while ^ is the scalar field from 
a particle in a circular orbit. The mass and charge of a 
particle are (i and q, and the mass of the black hole is m. 
The scalar field is "J"; however, we define ip = q^/n as a 
combination which occurs often in the description of the 
effects of a scalar field on the motion of a particle. In this 
perturbative analysis we always assume that <7 2 /V ^ m - 

In equations concerning the self force, expressions con- 
taining tp refer to the regularized field and must be eval- 
uated at the location of the particle. 



II. CONCEPTUAL FRAMEWORK 
A. Dissipative and conservative forces 

Mino [25] examines the self force on a particle in or- 
bit about a rotating black hole with a focus on how the 
Carter constant [26] evolves. Following Mino's ideas, 
one is naturally led to divide the self force into two 
parts which depend upon how each part changes under a 
change in the boundary conditions. One part is "dissipa- 
tive" and usually associated with radiation reaction. The 
other part is "conservative" . This distinction is most eas- 
ily described in terms of "Green's functions" which yield 
the parts of the field responsible for the dissipative and 
conservative parts of the self force. 

Let G ret be the retarded Green's function, which pro- 
vides the actual, physical field \E ,ret for the problem of 
interest, and let G be the Green's function [10] for 
the singular part of the field VE" . The regularized field 
vfr R = \]/ rct — vf/S provides the complete self force. The 
change in boundary conditions from outgoing radiation 
at infinity to incoming radiation is effected by using the 
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advanced Green's function G adv and its field i]/ adv , rather 
than the retarded quantities. 

The dissipative part of the self force changes sign under 
the interchange of the retarded and advanced Green's 
functions. Thus, the dissipative part of the field \I/ dls is 
uniquely determined by using a "Green's function" 

^ dls is a source-free solution of the field equation and is 
regular at the particle; no use of vf^ is required to find 
the dissipative part of the self force. 

The conservative part of the self force is invariant un- 
der the interchange of the retarded and advanced Green's 
functions. The conservative part of the field represents 
the half advanced plus half retarded field, but this is sin- 
gular at the particle and requires regularization by the 
removal of \P S . Thus, the conservative part of the regu- 
larized field is uniquely determined near the particle by 
using a "Green's function" 

G co„ = l G rct + l G adv_ G S^ (3) 

In a neighborhood of the particle, the resulting field is a 
source free solution of the field equation and is regular. 

A circular orbit about a Schwarzschild black hole pro- 
vides an example which clearly distinguishes the dissipa- 
tive from the conservative parts of the self force. For a 
particle in a circular orbit about a black hole, d t ^ and 
are both dissipative, as can be understood in terms 
of time-reversal invariance. We see below in Eqs. (16) 
and (17) that d t ^ and are responsible for remov- 
ing energy and angular momentum, respectively, from 
the particle at a rate which precisely balances the loss of 
energy and angular momentum outward through a dis- 
tant boundary and into the black hole. With t —> —t 
and ({> —> —cj) the motion of the particle is nearly identi- 
cal, only the boundary conditions on the scalar field are 
changed, and now energy and angular momentum enter- 
ing the system through the boundary are deposited on 
the particle. The dissipative self-force is small in this per- 
turbative analysis but its effect accumulates over time as 
the particle spirals slowly inward or outward, depending 
upon the boundary conditions. 

For a circular orbit d r ^ is conservative and, as shown 
below, provides a small addition to the centripetal force 
which affects the angular frequency via a change in the 
right hand side of Eq. (4). However, the angular fre- 
quency is unchanged under t —* —t and <fi — > — <j> which 
changes the direction of radiation flow at the boundary. 
d r ty and its effect on the angular frequency is indepen- 
dent of the direction of radiation imposed by the bound- 
ary conditions. With either boundary condition the con- 
servative self-force effect on the frequency is small but the 
effect on the phase of the orbit accumulates over time. 

If a particle is not in a circular orbit, then these sim- 
ple relationships between the components of d a ij and 
the conservative and dissipative parts of the self force 
no longer hold. 



B. Stability of a circular orbit 

The notion of the stability of a circular orbit in the 
context of the self-force warrants further discussion. 

A stability problem presupposes the existence of a me- 
chanical system in equilibrium. A perturbation analysis 
of the system's small oscillations often reveals that the 
natural frequencies of the system are complex eigenval- 
ues dependent upon some set of parameters. The sign 
of the imaginary part of a natural frequency determines 
whether the amplitude of a small oscillation grows, di- 
minishes or stays constant. The simplest cases are when 
the frequency is either purely real or purely imaginary. 
An imaginary frequency with an appropriate sign then 
signifies an unstable mode. 

For the system of a small particle of mass fi orbiting a 
larger black hole of mass m, the particle moves along a 
geodesic, if /i is considered infinitesimal. In this case it 
is well known that for a circular orbit at Schwarzschild 
coordinate r = r Q the angular frequency fi , with respect 
to Schwarzschild coordinate time, is given by 



If the orbit has a small eccentricity then elementary anal- 
ysis [27] reveals that the frequency fl r of the radial oscil- 
lations is determined by 

n?= ™(r„-6m). (5) 

If r Q is very large, then f2 Q and £l r are nearly equal, and 
the slightly eccentric orbit is an ellipse obeying Kepler's 
laws. The difference between Q and fi r leads to the rate 
of the precession of the ellipse, 

and is directly responsible for the relativistic contribution 
to the precession of the perihelion of Mercury. For an 
orbit just inside 6m, the two solutions for Q r in Eq. (5) 
are both imaginary and one corresponds to an unstable 
mode. The orbit at r Q = 6m is the innermost stable 
circular orbit — the ISCO. 

For realistic boundary conditions, the expectation is 
that at a large separation the emission of gravitational 
radiation from a binary circularizes the orbit [28] and 
causes the two objects to spiral slowly towards each other. 

When r Q is outside but comparable to 6m, gravita- 
tional radiation evolves the orbit slowly inward on a 
secular timescale m/fj, times longer than the dynamical 
timescale, with the rate of inspiral dr a /dt — 0(n/m). 
The particle makes a transition to a plunge [29] when 
it reaches the ISCO, and the plunge occurs over a dy- 
namical timescale. The resulting gravitational wave form 
appears as a sinusoid with a slowly increasing frequency 
until the angular frequency of the ISCO is reached. Then, 
after a brief burst from the plunge, the wave form is de- 
termined by the frequencies of the most weakly damped 
free oscillations of the black hole [30-32] . 
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The secular evolution caused by gravitational radiation 
keeps the particle from ever being in true equilibrium, 
and the question of the stability of a circular orbit may 
seem ill posed. Nevertheless, the gravitational wave form 
changes character at a frequency near that in Eq. (4) for 
r = 6m even though the particle is never actually in 
equilibrium. The dependence of the transition frequency 
upon the scalar field self force is a major focus of this 
manuscript. 

To form a well-posed problem related to the ISCO and 
the angular frequency of the transition to a plunge, we 
consider an unphysical system with boundary conditions 
which have equal amounts of outgoing and incoming radi- 
ation. Such a system must be constructed by use of a reg- 
ularized half-advanced and half-retarded Green's func- 
tion. In this case, the dissipative force vanishes entirely 
while the conservative force and its effects upon f2 Q and 
Q r are unchanged from the case with realistic outgoing- 
wave boundary conditions. With equal amounts of out- 
going and incoming radiation, equilibrium configurations 
exist, and stability analyses may proceed. 

In Section V we show that the conservative force from 
V a ^ changes the right hand side of Eq. (5), and therefore 
the location of the actual ISCO by a fractional amount of 
order q 2 / [im for a scalar field; for gravity this fractional 
amount would be of order fi/m. 



III. SCALAR FIELD EFFECTS ON A 
PARTICLE'S MOTION 



The normalization of u a implies that 



vTu a 



(E 2 - R 2 ) .P_ 
1 - 2m/ R + R 2 



or 



E z - R 2 = (1 - 2m/R){\ + J 2 /R 2 ). 



(11) 



(12) 



B. Scalar field modifications of the geodesic 
equation 

Quinn [33] considers the interaction of a scalar field 
^ with a particle of constant bare mass /j,q and constant 
scalar charge q. He carefully demonstrates that the equa- 
tion of motion is 

u b V b (fiu a ) = (iu a + fiu b V b u a = gV a *, (13) 

where [i = fio — q*&- This result can be obtained from a 
general action principle. In what follows, is assumed to 
be the regular field 5» R evaluated at the particle. Thus, 
the quantity /U changes when the value of the scalar field 
at the particle changes. The fractional change in /j, is 



(14) 



which is shown in Figs. 1 and 2 for circular orbits of the 
Schwarzschild geometry. 

The projection of Eq. (13) orthogonal to u a yields the 
acceleration from the self-force 



A. Description of motion 

The functions (T(s), R(s), tt/2, <&(s)) describe a parti- 
cle's worldline in the equatorial plane of a black hole in 
Schwarzschild coordinates (t,r,0,(j>). The proper time s 
is measured along the worldline, and the four velocity 



u a = dx a /ds 



(7) 



is normalized to unity, u a u a = — 1. Additionally, the 
worldline and the Killing vector fields t a d/dx a = d/dt 
and (j) a d/dx a = d/d(j) define 



E 



-t a u a and 



(8) 



which, for geodesic motion, are the energy and angular 
momentum, per unit rest mass, respectively. The com- 
ponents of u a are thus 



E 



1 - 2m/ R 



R, 0, J/R 2 , 



(9) 



where the dot denotes a derivative with respect to s. 
For any worldline the angular frequency, with respect to 
Schwarzschild coordinate time, is 



« = S/T 



d$ J(R - 2m) 



dT 



ER 3 



(10) 



u b V b u a = (q/ti)(g a b + u a u b )V b ->S 



(15) 



which modifies the worldline of the particle through 
spacetime. In our perturbative analysis, we only con- 
sider the first order effects of the scalar field. The right 
hand side of Eq. (15) is explicitly first order in ^, and 
the change in /x has only a second order effect on the 
acceleration. Thus, for the purposes of describing the 
world line of the particle we treat /i as constant, and it 
is convenient to define ip = q^ / /J-. For simplicity we also 
assume that the scalar field is symmetric under reflec- 
tion through the equatorial plane and that the particle 
remains in the equatorial plane. 
The components of Eq. (15) yield 



dE/ds = E= -d t ip + Eu b V b ip 



dJ/ds = j = d^ip + Ju b V b ip 



(16) 



(17) 



and 



R 



m{E 2 -R 2 ) R-2m 2 



R(R - 2m 
R-2m 



R 



R 4 

d r i> + Ru b \7 b ip. 



(18) 



The normalization (12) is a first integral of the equa- 
tion of motion, and Eq. (18) follows directly from Eqs. 
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FIG. 1: The regularized field ^ R m/q at the particle and the 
radial component of the self force (<9 r $ R )m 2 / q as a function 
of the radius r Q for circular orbits close to the black hole. 



(12), (16) and (17), which form a complete set of equa- 
tions describing the equatorial motion of a particle inter- 
acting with a scalar field. 

Together, Eqs. (12) and (18) imply that 



R ~ 'I? 



R — 3m 
R 4 



J 



+ — d r ip + Ru b V b iP, 



(19) 



which is convenient for analyzing slightly eccentric orbits. 

For a particle in circular motion, d r ip is independent of 
the direction of radiation at the boundaries and provides 
a conservative force. However, d t ip and d^ip on the right 
hand sides of Eqs. (16) and (17) are dissipative compo- 
nents of the self force and change sign if the direction of 
radiation at the boundary is changed. In fact, the conser- 
vation of energy and angular momentum guarantee that 
Eqs. (16) and (17) are consistent with the flow of energy 
and angular momentum across the boundaries. This is 
manifest by matching appropriate independent solutions 
of the source-free wave equation at the orbit of the par- 
ticle. If a particle is not in a circular orbit, then these 
simple relationships between the components of d a ip and 
the conservative and dissipative parts of the self force no 
longer hold. 



IV. SELF FORCE EFFECTS ON CIRCULAR 
ORBITS 

The scalar field ip of a charged particle in a circu- 
lar orbit rotates with the particle and has a symmetry 
described by 



d t ip + Q d<j,ipo = 0. 



(20) 



In particular, this implies that w a V a ?/>o = at the par- 
ticle and simplifies the description of the self force for 




FIG. 2: The same as Fig. 1 except for larger radii. 



circular motion. A worldline is "instantaneously circu- 
lar" at radius r Q if R = R = 0. In this case 



J 



m 

3m \ ^ 2 



2m 



r Q — 6m \r. 
follows from Eq. (19), and 



E 2 = 



(To - 2m) 2 
r (r - 3m) 



d r 1p 



(1 - r d r ip ) 



(21) 



(22) 



then follows from Eq. (12). The angular frequency is 
given by Eq. (10) 



2 171 r o — 3m 

° ~~ 3 2 



cWo + 0(V> 2 ). 



(23) 



These reduce to the usual expressions for the circular 
geodesies of the Schwarzschild geometry when the scalar 
field ip a is removed. The resulting fractional changes in 
J, E and f2 Q caused by the self force for an orbit at radius 
r n are 



AJ 



r (r - 2m) 
2m 



d r tp , 



and 



A0 o _ r (r - 3m) 
il 2m 



d r ip . 



(24) 



(25) 



(26) 



Rcfs. [19, 23, 34] give detailed descriptions of the calcu- 
lation of the retarded scalar field for a particle in a circu- 
lar orbit of the Schwarzschild geometry. Rcfs. [12, 15, 17- 
21] describe the evaluation of the required regularization 
parameters for the self force. In Appendix B, we give the 
regularization parameters for the scalar field from 
Ref. [20]. Refs. [19, 23] describe implementations of the 
Barack and Ori regularization procedure for the numer- 
ical evaluation of the self force for circular orbits of the 
Schwarzschild geometry. 
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30 


-3.661710e-5 


7.171924e-8 


-1.54904e-8 


5.08927e-6 


-3.82435e-4 


50 


-7.889518e-6 


6.346791e-9 


-8.53532e-10 


6.44498e-7 


-1.07242e-4 


70 


-2.877222e-6 


1.284529e-9 


-1.26087e-10 


1.66932e-7 


-4.65496e-5 


100 


-9.884245e-7 


2.356504e-10 


-1.65135e-ll 


4.00531e-8 


-1.92277e-5 


200 


-1.239866e-7 


8.642538e-12 


-3.11790e-13 


2.51466e-9 


-3.44412e-6 



TABLE I: A selection of our computed values of the regularized scalar field evaluated at the particle and the radial component 
of the self force, for circular orbits of radius r . Also shown are F r , Ft and F<p, defined in Eqs. (33)-(35), for slightly eccentric 
orbits with r D > 6m; for r D = 6m, their limiting values are given; for r D < 6m slightly eccentric orbits do not exist and F r , F t 
and F$ are not defined. 
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FIG. 3: The fractional changes in J, E and fi , from the self 
force as a function of the radius r for circular orbits close 
to the black hole. These quantities are calculated using Eqs. 
(24)-(26). 
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Figs. 1 and 2 show the regularized ip R and the radial 
component of the self force d r ip R evaluated at the particle 
as functions of the radius of the orbit. Figs. 3 and 4 show 
the fractional changes in J, E and tt from the radial 
component of the self force. Table 1 gives a selection of 
our numerical calculations of both ip R and d r ip R at the 
particle as a function of R. The force is outward, falls 
off approximately as r~ 5 , as noted by Burko[23], and if> R 
falls off as r Q ~ 3 for large r . Each of Cl , E and J for a 
circular orbit is decreased by the self force. 



V. SELF FORCE EFFECTS ON SLIGHTLY 
ECCENTRIC ORBITS 

A particle in an orbit with small eccentricity, in the 
Schwarzschild geometry, has periodic motion in the ra- 
dial direction with a frequency fl r given in Eq. (5) . In this 
section we describe the effects of the self force on Q r and 
upon the rate of precession of the orbit f2 pr . We require 
that the scalar field ip have boundary conditions with 
equal amplitudes of incoming and outgoing radiation at 
both the event horizon and infinity. This simplifies the 
discussion of fi r and allows us to pose a well defined sta- 
bility problem for the ISCO. The conservative "Green's 
function" G con determines the regularized field. 



A. Slightly eccentric geodesies 

In this perturbative analysis, a slightly eccentric 
geodesic about a fixed radius r with a small amplitude, 
SR < q 2 jix < r a q 2 / {i 2 < r Q , is described by 



FIG. 4: The same as Fig. 3 except for larger radii. 



R{T{s)) =r + SRcos(n r T{s)) (27) 
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and 

$(T) = Q T+^^sin(fyT). (28) 
ait s i r 

The remainder of the analysis treats SR as a small quan- 
tity, and only terms through first order in SR are retained. 
The angular frequency of this geodesic orbit is a function 
of time, 

fy = ^i = + ^fficos(O r T). (29) 

The frequencies f2 and O r are given in Eqs. (4) and (5). 
The quantity dQ^/dR represents the change in angular 
frequency, from Eq. (10), with respect to a change in 
radius while E and J are held constant, 



dn^ _ 2{R-3m) 
~~dR ~ ~ R(R-2m) 

The radial velocity is 



0,4. 



d Ft 

— = - fl r 6Rsm(fl r T). 



(30) 



(31) 



B. Self force effect upon J for slightly eccentric 
orbits 

For slightly eccentric orbits, with self force effects in- 
cluded, J is not a constant of the motion, even with no 
dissipation. The effect of the self force on J is given in 
Eq. (17), 

dJ/ds = d<pip + Ju a S7 a ip or 
dJ/dT = (u^d^j + Jidt + fl^iP 



(32) 



The ip for a slightly eccentric orbit is described in Ap- 
pendix A. We note here that it has two parts. The larger 
part ip is equal to the field which would result from pure 
circular motion at r Q and consists of frequencies which 
are integral multiples of Q - For a circular orbit, in gen- 
eral, (dt + ^l d^) ip = 0, and with conservative boundary 
conditions d t ip = d^ipo = at the particle. The smaller 
part is proportional to SR. In Appendix A we show that 
for slightly eccentric orbits an expansion for ip around a 
circular orbit at r a gives 

d t tp = -F t Q r SRsin(Q r T), (33) 



and 



dfip = -i^fi r ffisin(fi r T), (34) 



d r i\) = [d r ip ] ro + F r 5Rcos(n r T), (35) 



where F t , F^ and F r depend only upon r and, are inde- 
pendent of both SR and t. In this section the subscript 



r Q on [d r ip ]r implies that d r ip is to be evaluated at 
the circular orbit r and not at the actual location of the 
particle. In Eq. (32) the coefficient of d r ip is already first 
order in SR, and it is sufficient to use only the circular 
orbit value [d r ip ] ro at the particle. Thus, 



dJ dJ fdR*- 
dR ~ dT \dT / 

= JF t + EF4 + J[dri) \ ro , 



(36) 



where we have used Eqs. (31), (33) and (34), along with 
Eqs. (9)-(12). 



C. Self force effects upon Q. r for slightly eccentric 
orbits 



The exact radial equation of motion (19), which in- 
cludes the self-force, is 

» m R— 3m 2 

R = ~& + — u^ J 

+ R ^ m d r ^ + Ru b V b 4>. (37) 
R 

For describing slightly eccentric orbits, we expand this 
equation around the circular orbit at r a by letting R — > 
r Q + SRcos(fl r T) and dropping all terms of order SR 2 . 
The 0(SR ) part is equivalent to Eq. (21) for the circular 
orbit value of J 2 . The Ru b V b ip term is 0{SR 2 ). The part 
of Eq. (37) which is first order in SR is 



R = 



dR 



to R — 3m 



+ 



R 2 R 4 

R-2m " 
-o r ip 



J 



R 



5Rcos(n r T). (38) 



After use of Eq. (35), which gives d(d r ip)/dR = F r , this 
becomes 



R 



m(r -6m) 3(r Q - 4m)(r - 2m) 

+ — [d r 1po\r a 



r3(r - 3m) 



r 2 (r Q - 3m) 



2m (r G - 3m) dJ 2 



| (r Q - 2m) 



Eq. (27) implies that 



R=- 



r4 dR 
SRcos(fl r T). 

Q 2 .E 2 5Rcos(n r T) 
(1 - 2m/r Q ) 2 



(39) 



(40) 



Using Eq. (22) for E 2 , we finally obtain 



m 2 J 
= -^(r -6m) E (r -3m) 2 (JF t + EF <j> ) 



-(r - 2m)(r G - 3m){3[d r ^ } ro + rF r )(41) 
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FIG. 5: For q 2 /fi <C r a — 6m, the change in fi r from the self 
force as a function of the radius r for slightly eccentric orbits 
close to the black hole. The limiting behavior at q 2 /fi -C 
T — 6m « m is given in Eq. (46). See Eq. (43) for g 2 //i ~ 
r — 6m. 



through first order in i[>. This provides the scalar- field 
self- force correction to Eq. (5) for f2 2 . 

To examine the effect of the self force on Q r for r w 
6m, let / represent all but the first term on the right 
hand side of Eq. (41). The change in O 2 caused by the 
self force is then 



FIG. 6: From numerical analysis, the fractional change in £l r 
from the self force as a function of the radius r for slightly 
eccentric orbits far from the black hole. 



The right hand sides of Eqs. (43) and (46) are equal at 

£--l«J^/ = 0.00511255 q 2 /,m. (47) 

At larger separations r — 6m > m, AQ r is still approx- 
imated by Eq. (45) and scales as l/r Q 2 as illustrated in 
Fig. (6). ' 



A(Q 2 ) = / = / + /i(r o -6m) + 



(42) 



Numerical analysis determines the values /o = 9.46768 x 
ICT'Y/mw 3 , and fx = -3.2318 x 10~ 5 q 2 / fj,m 2 . For r G = 
6m, 



AQ r = fl /2 = 9.73020 x 10" 3 g/ y /2 m 3/2 . (43) 

f -»6m 

More generally, the change in Q r caused by the self 
force is 



rn 



—r{r - 6m) + / 



1/2 



-^(r -6m) 

' O 



1/2 

(44) 



When r Q increases away from 6m, A£l r decreases and 
changes scale 

fr 2 

for r Q - 6m > q 2 /n, Ail r w J —^- [m(r G - 6m)]~ 1/2 



= 0(q 2 /iim 2 



(45) 



The numerical calculation of AQ r , based upon Eq. (45), 
is presented in Figs. (5) and (6). Fig. (5) illustrates the 
small end of this range, where q 2 /fi <C r a — 6m <C m and 
/ w /o- In this case, 

AO. w 6.95730 x 10" 4 (r o /6m- iy 1/2 q 2 / (im 2 {46) 

gives the limit of the curve in Fig. (5) as r a approaches 
6m. 



D. Self force effects upon fi pr for slightly eccentric 
orbits 

The rate of precession of a slightly eccentric orbit 



(48) 



is not particularly elegant when written in terms of m, 
r Q and the components of the self force. However, 



rn 



1/2 



(6m) 3 / 2 



/ 1/2 for r = 6m. (49) 



In general, 



Afi pr = An o - An r 



(50) 



where A0 o and AO r may be obtained from Eqs. (26) 
and (44) . 

Fig. 7 illustrates AQ pr for orbits for smaller values of 
r Q , and the curve for Af2 pr has the same limit to that of 
— Ail r in Fig. (5). The fractional change Af2 pr /f2 pr from 
the self force is shown in Fig. 8 for larger values of r Q 
where Af2 pr /f2 pr scales as r Q _1 . 

We summarize the falloff at large r for a variety of 
quantities in Table II. 

Our earlier manuscript [19] described the numerical 
evaluation of d r tp^ for a circular orbit in great detail. 
The analytically known regularization parameters (the 
multipole moments of d r ip s ) were subtracted from the 
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FIG. 7: For q/m «r„- 6m, the change in f2 pr , from the self 
force as a function of the radius r a for slightly eccentric orbits 
close to the black hole. 



numerically determined d r tp^ m . A few additional reg- 
ularization parameters were determined numerically and 
also subtracted from d r tp^ m . The remainder was summed 
over I up to about 40. 

The main difficulty revolved around the evaluation of 
drtfj^ 71 with sufficient accuracy that the final sum gave us 
good precision. Starting with approximately 13 signifi- 
cant digits for d r ipo, after the regularization parameters 
were subtracted about 8 significant digits remained. The 
main numerical task for evaluating the scalar field self- 
force effects on fi r and Q pr is very similar to this earlier 
work. A significant difference, however, is the need to 
compute F t , and F r , introduced in Eqs. (33)-(35). 
The details required for determining the part of the field 
which depends upon the slight eccentricity of the orbit 
are described in Appendix A, and the regularization pa- 
rameters are given in Appendix B. 



VI. SELF FORCE EFFECT ON THE ISCO 

The innermost stable circular orbit is characterized as 
that orbit for which £l r is zero. The self force changes 
the radius of the ISCO from 6m by AR where 

AR is = (i? is -6m) (51) 
= 180m 2 d r Vo + {3/2)dJ 2 /dR + 432m 3 F r . 

This result follows from equating the coefficient in square 
brackets in Eq. (39) to zero, replacing r by 6m in all 
terms which are 0(q 2 /[im) and solving for (r Q — 6m). The 
self-force correction to the angular frequency of the ISCO 
is given by Eq. (23) when evaluated at r Q = 6m + Ai? is , 

Clf s = r ~ 4 m(r - 3Ai? is ) - r ~ 2 (r - 3m)d r ip 
= m(6m)~ 3 - 3mAi? is (6m)~ 4 

- (12m)- 1 d r ip . (52) 



0.05 




10 ~ 100 

r /m 



FIG. 8: For q/m <ti r — 6m, the fractional change in Q pl 
from the self force as a function of the radius r a for slightly 
eccentric orbits far from the black hole. 

TABLE II: Falloff at large r D for a variety of interesting quan- 
tities involving the scalar-field self-force. 



Quantity 


behavior 




~ r ~ 3 


drip* 


~ r ~ 5 

~ / o 


flo 


~ r ~ 3 / 2 

~ / O 




~ / o 




' o 


AQ, r 


~ r ~ 7 / 2 

1 O 




~ r ~ 5 / 2 

~ / o 


Af2 pr 




E 


~ r o 


AE 




J 


~ ro 1 / 2 


AJ 


~ r ~ 5/2 



The fractional change in fi; s from the scalar field self- 
force is 

AQ, is /fl is = -Ai? is /4m - 9md r i/j . (53) 

A summary of our numerical results for the effect of 
the scalar field self-force on the ISCO is given in Table 
III. Our numerical work primarily followed that for eval- 
uating O pr . However, one subtlety involved the need for 
evaluating both F t and F^ in the limit that r — ► 6m. 
For that task we found each quantity at approximately 
thirty points between r = 6m and 7m. A polynomial 
fit to these data, using a variety of subsets of the thirty 
points and different numbers of polynomial coefficients, 
provided robust results for the required limits. 

VII. DISCUSSION 

Two parts of our analysis of the self force effects on 
slightly eccentric orbits are algebraically taxing while not 
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TABLE III: Quantities of interest regarding the ISCO. 



AR is x [i/q 2 


-0.122701 


AQi B /Qi B x fim/q 2 


0.0291657 



particularly difficult conceptually. These are the calcula- 
tions of the regularization parameters and the matching 
of the homogeneous solutions of the field equation across 
the orbit of the particle. While these two steps are in- 
dividually challenging, our confidence in the ultimate re- 
sults is bolstered by the appropriate convergence of the 
sum over the Im modes of the components of the self 
force. An error in cither the analytical or the numerical 
work involving either the matching or the regularization 
parameters would be immediately heralded by a lack of 
convergence of the sum over modes. Consequently, we 
deem these results trustworthy. 

We see that the effect of the self force from a 
scalar field on the innermost stable circular orbit of the 
Schwarzschild geometry is to move the ISCO inward and 
to increase its angular frequency. It is tempting to gen- 
eralize this result to the gravitational self force case. But 
we will not do so. There is a significant difference be- 
tween the self force effects of a scalar and a gravitational 
field. In particular, as we mentioned in the introduction, 
for gravity the self force already has an important effect 
at the Newtonian level [4]. Namely it is responsible for 
ensuring that the particle and the black hole both or- 
bit their common center of mass. This is not the case 
for the scalar field. The difference may be traced back 
to the fact that a black hole has no scalar charge. The 
charged particle responds only to the gravitational inter- 
action with the black hole and the scalar field interaction 
with its own ^> R , and the particle motion deviates from a 
geodesic of the Schwarzschild geometry only because of 
its scalar field. In this limit, the black hole is not affected 
at all and remains fixed in space. And the particle orbit 
is centered upon the black hole, not upon the common 
center of mass. To see the motion about the center of 
mass, it is necessary to consider the gravitational self- 
force problem which we will return to in a later paper in 
this series. 

A well defined formulation of the stability of the ISCO 
requires the imposition of equal amounts of outgoing and 
ingoing radiation at the boundaries to make the system 
dissipation-free. However, actual slow inspiral into a 
black hole has only outgoing radiation, and one might 
wonder about the relevance of our calculation to the ac- 
tual physical system. Ori and Thorne [29] outline a care- 
ful treatment of the actual slow evolution of a small ob- 
ject inspiralling and making the transition to a plunge 
into a black hole. They also indicate how self force ef- 
fects could be included in their analysis. An extension 
of our results in the manner outlined by them is beyond 
the scope of this manuscript, but may be returned to in 
the future. 
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APPENDIX A: SOURCE DECOMPOSITION FOR 
NEARLY CIRCULAR ORBITS 

In this Appendix m refers to the spherical harmonic 
index while M refers to the mass of the black hole. Else- 
where in this manuscript m is used for both of these 
quantities without confusion. 

The effects of the self force on il r and on the ISCO 
are governed by the scalar field tp for a slightly eccentric 
geodesic. Finding the field in this case is more difficult 
than for a circular orbit. Apostolatos, et al.[28], consider 
the stress-energy for a point mass moving along such a 
geodesic, in their Sections IV. B and IV. C, and provide 
an expansion of the source in powers of a small constant 
SR. Their analysis is easily modified for a scalar charge. 

A slightly eccentric geodesic in the equatorial plane is 
described by 

R(t) =r + 5Rcos{Cl r t) (Al) 

and 

$(t) = n a t + s$(t) = n t + c ^±^ S in(n r t), (A2) 

dR il r 

where the angular frequency of the orbit is 

f^ = ^=f2 + ^fficos(Q r t). (A3) 
dt dR 

The frequencies Q and Sl r are given in Eqs. (4) and (5). 
The quantity dVL^/dR represents the change in angular 
frequency with respect to a change in radius 

d% 2(r Q -3M) o 

llR ~ -r (r -2M) n *' (A4) 
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while the conserved quantities E and J are held constant an integration by parts over the angle <f), yields 
at their circular orbit values, 



E 2 = (r -2mY 



J 2 = 



(A5) 



f (v — 3m) r Q — 3m 

The scalar-field source which moves on such an orbit 



Qe m (t,r) = <j> gYl 



q(r a 2M ) y/m( ^ A0)e -^V 



is 



rjE 
1 — im 



(-gyWSlr-RiaMt-Zia)] 



sm(il r t) 



dR n r 

x [6{r - r Q ) - 5Rcos(Sl r t)5'(r - r Q )] . 



(A7) 



x 5(6 - tt/2) S[t - T(s)] ds 
= 5[r-r -SRcos(n r t)] 



x 5 



dQ.^ SR . ' 

~~dRrT sm ^ ' 



after an integration over s and substitutions from Eq. 
(Al) and (A2). An expansion of the delta functions 
about the orbit, with SR being small, along with a spheri- 
cal harmonic decomposition of this source, which includes 



This reveals that the source has a frequency spec- 
trum consisting of the harmonics of the angular fre- 
quency io m = mfl along with sidebands at frequencies 
, /9 s ! (\a\ w m = UJm ^ ^r- The amplitude of the sidebands are 
^ 77 ' ^ ' proportional to SR. 

The scalar field wave equation is 



V 2 * = -4tt£>. 



(A8) 



The separation of variables of VE" yields 



(A9) 



im 



imuj 



where I = . . . oo, m = —£...£ and u) — {uj m , u) m , The radial equation for the r dependence of any I mode with 
a frequency u is 



d 2 *? m (r) , 2(r-M) d*? m (r) 



c?r 2 



+ 



r(r - 2M) dr 



+ 



w 2 r 2 



(r - 2M) 2 r(r - 2M) 



1 - 2M/r 



(A10) 



For self force calculations it is convenient to use tp = (q/fj,)^ and to divide each ipg m into the u m part ip% m and the 



side band parts 8R\ 



im' 



(All) 



The numerical determination of these parts of ip£ m re- 
quires the source-free solutions of Eq. (A10) with appro- 
priate boundary conditions, and then the proper match 
of these solutions across the orbit of the particle at r . 
The matching conditions, from Eqs. (A7) and (A10), are 



Wtm\ ro = [^ m (r + e)-V? m (r -e)] 



dr h 



tw 



dr 



(A12) 
(A13) 



dxt 



dr 



Anq 2 
~Hr Q 2 E 



y; m (7r/2,o) 



M 



± 



m dCl d 



r {r - 2M) 2Q, r dR 



(A14) 



(A15) 



where, in the Appendix only, [ ] ro on the left hand side 
denotes the discontinuous change in a quantity across the 
orbit at r . 



The scalar field tp is a purely real field, and it is convenient to combine the m and — m contributions 

^tm + 1>t,-m = Wt m + VC) cos[m(0 - n o t)]Y em {9, 0) + ( X | m + Xi^SR cos[m(<£ - n o t) - n r t]Y em (6, 0) 
+ {X7 m + XlJSRcosim^ - Slot) + n r t]Y em (6, 0), 



(A16) 
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where we have used the fact that ip is the conservative field given by ^(ip rct + i/> adv ). The contribution from m = 
for a given I is just one half of that in Eq. (A16) with m set to zero, 



(A17) 



Note that xJ,q = xto f° r the conservative field. 

In Section VI we require ^ m + tpe,- m and its derivatives evaluated at the particle. The location of the particle is 
given in Eq. (Al) and (A2). For m ^ 0, an expansion about r , retaining only terms through first order in SR, yields 



{ipem + ipe-m). 



1 + 



x m m + «J cos[m(0 - n o t)]Y tm ) 



+ (Xtm + Xlm + Xi m + X im ) SR COs{Sl r t)Y erl 



(A18) 



where the subscript p implies evaluation at the particle. From here through the remainder of this Appendix A, the 
spherical harmonic Yg m is to be evaluated at (ir/2,0), and the other terms on the right hand sides are to be evaluated 
at (r, <j>) — (r a , il t) only after all appropriate derivatives have been taken. After simplification, this becomes 

(1> tm + ^,-m) p = (V> £ ° m + ^)Y tm + {d r r tm + dr^Z + xL + xZ + Xj m + X7*) SR cos{Q r t)Y lm . (A19) 

Similar expansions start with Eq. (A16) and ultimately provide the 4> derivative 

d 



(ll>tm + C) + m(Xi m + Xlm - X £m ~ Xtot) 



SRsm(n r t)Y tm , (A20) 



the t derivative 
d 



£l r dR 



(V>?m + ViL) ~ W+(xL +Xtm)+ "miXem + XeZ) 



5Rsm(n r t)Y em , (A21) 



and the r derivative 



+ d rX~lm + d rXtm + d rXj m + d r Xj£)SRcOs{Vt r t)Y lm . 



(A22) 



In Section VI we require the sum over I and m of these 
three previous derivatives. Accordingly, we define F t , 
and F r from 

(dti> R ) p = -F t n r 5Rsm(n r t) (A23) 
(d V R ) P = -F^arSRsmiQrt) (A24) 

and 

(d r ^ K )p = (d r ipo)r + F r 5Rcos(Q r t), (A25) 

where the required regularization is described in Ap- 
pendix B. 

APPENDIX B: REGULARIZATION 
PARAMETERS FOR tp 

We describe the regularization of the scalar field as de- 
veloped by Barack and Ori [12, 17] for a particle in a 



circular orbit of the Schwarzschild geometry. Our nota- 
tion follows that of Refs. [19-21] . All of these are required 
reading for a thorough understanding of this Appendix. 

The scalar field is regularized at the location of the 
particle by subtracting the singular part of the field tp s 
from the actual field tp. The remainder -0 R = -0 — -0 s 
is then guaranteed [10] to be a regular solution of the 
vacuum scalar field equation in the vicinity of the parti- 
cle, and the derivatives of tp K at the particle provide the 
required self-force resulting from the particle interacting 
with its own field. 

The mode sum regularization procedure [12, 17] de- 
scribes the multipole decomposition of ip R in terms of 
the decompositions of tp and s , 

= J2^™(t,r)Y em (6,t) 

Im 

= bim{t,r) - VL(*,r)] Y tm (6,<f>). (Bl) 

The numerical determination of ip trait, r ) for a slightly 
eccentric orbit is discussed in Appendix A and in 
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Ref. [19]. ip s , however, is singular at the location of the 
particle, and only well defined in a neighborhood of the 
particle. Nevertheless, its multipolc decomposition over 
a two-sphere of radius r is finite, even if r coincides with 
the the radial coordinate of the particle. The decompo- 
sition is not unique because of the ambiguity in the defi- 
nition of ip s away from the particle. However, the mode 
sum regularization procedure remains valid because its 
sum is only required in a neighborhood of the particle, 
where the sum must equal ip . Thus, in evaluating ■0 R 
and its derivatives at the particle the individual ipe m m 
Eq. (Bl) are not unique but the sum in Eq. (Bl) and its 
derivatives converge to unique values. 

Barack and Ori [12, 17] find it convenient to do the 
sum over m first and then to describe the multipole de- 
composition of a derivative of ^ S as 



+ o{r 2 ) 



(B2) 

where the 0(£ 2 ) terms yield precisely zero when 
summed over £. The constants A a , B a and C a arc in- 
dependent of £, and are determined by a multipole de- 
composition of an expansion of -0 s about the location of 
the particle. 

The required regularization parameters for the deriva- 
tives of dttp are derived from Eqs. (8a)-(8d) of Ref. [17]. 
We discuss only those parameters which have not previ- 
ously appeared in an actual application[19, 23]. In our 
notation, with R representing a derivative of R with re- 
spect to proper time s, these regularization parameters 



are 



A ±t = ± 



A± r = =F 



i 2 R 



p(R 2 + J 2 ) 



q 2 E 



p,R 2 (E 2 - R 2 ) 



A ±(j> = 

q 2 ERR , > 
Bt ~ 2p{R 2 + J 2 f/ 2 {Fl/2 ~ 2F - 1/2) 



(B3) 

(B4) 
(B5) 
(B6) 



q 2 R 2 



(2E 2 - R 2 )F 1/2 - (E 2 + R 2 )F_ 1/2 



2p(R-2M)(R 2 + J 2 ) 3 / 2 

= q 2 RR(F 1/2 - F_ 1/2 ) 
2pJ(R 2 + J 2 ) 1 / 2 



Ct — C r — Ca 



0. 



(B7) 
(B8) 

(B9) 



Here the hypergeometric function is represented by F q 
2 Fi(q, 2', 1; z) where the argument z = M/(R — 2M). 



These parameters may be expanded by use of Eq. (Al), 
which implies that 



R = - 



5REtt r sm{fl r t) 
1 - 2M/R 



(BIO) 



Through first order in SR the non-zero regularization pa- 
rameters for a slightly eccentric orbit are 



A± t = T^—^SRsm(n r t) 

/j,r 2 



lxr 2 E 



1 SRcos(n r t) 



(Bll) 



(B12) 



B r — — 



q 2 n r (Fj 



1/2 



2F_ 



1/2 



2/zr (r 2 + J 2 )i/2 



5Rsin(n r t) (B13) 



q 2 (2F 1/2 



F_ 



1/2 



) 



2fir (r G 2 + J 2 ) 1/2 
2r Q 2 + J 2 



1 - 



r (r 2 + J 2 ) 



SR cos(51 r i) 



q 2 M (2F[ 



/2 



F' 



1/2 



2^r (r - 2M) 2 (r Q 2 + J 2 ) 1 ^ 



SR cos(O r i) 



(B14) 



Ba = — 



q 2 n r (F 1/2 - F_ 1/2 ) 
2fiJ(l - 2M/r o y/ 2 



SRsin(n r t) (B15) 



In this expansion, F[, 2 is the derivative of the hyper- 
geometric function F X j 2 with respect to its argument z. 
Both F 1/2 and F[, 2 are evaluated at z = M/(r a - 2M). 

The regularization parameters for F t and Fa,, defined in 
Eqs. (A23) and (A24), are obtained by removing the fac- 
tor — Q, r SRsm(Q, r i) from A± t , B t and Ba,. Similarly, the 
regularization parameters for F r , defined in Eq. (A25), 
are obtained by removing the factor SR cos(Q r t) from the 
SR terms of A± r and B r . 

The regularization parameters for the scalar field, 
alone, warrants further discussion [20]. In a particular 
locally-inertial coordinate system (T, X, Y, Z) , the singu- 
lar field near a scalar charged particle is simply 



lP S = q/p + 0(p 3 /ll 4 ) 



(B16) 



where 1Z is a length scale of the geometry in the vicin- 
ity of the particle, and p 2 = X 2 + Y 2 + Z 2 . For 
the special case that the particle is in a circular orbit 
about a Schwarzschild black hole, a coordinate transfor- 
mation between the special (T, X, Y, Z) coordinates and 
the usual Schwarzschild coordinates allows p to be writ- 
ten as a function of Schwarzschild coordinates, and the 
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expansion of 1/p about p = is given in Eq. (6.22) of 
Messaritaki [20] . The terms of interest are 



1 

9 



-+6 1 



3m 



8rl(r Q - 2m) \x 



(r + m) 1 
r a X 2 



(B17) 

where . . . refers to terms which vanish as r — ► in such 
a manner that they have no contribution to the regular- 
ization parameters. In this equation we use 



9 2 



where 



r (r - r Q f 
r Q — 2m 



2r 



2m 



° r Q — 3m 
m sin 2 $ 



x(l-cos6) (B18) 



2m' 



(B19) 



The angles (0, $) are derived from a rotation of the 
Schwarzschild coordinates which puts the particle on the 
9 = axis [19]. 

We find that, through the e 1 term of Eq. (B17) 



Vlofr = r Q ) = B n 



2V2LV 



{21- l)(2£ + 3) 



(B20) 



where and result from the e _1 and e 1 terms of 
Eq. (B17) respectively. Appendices C and D of Ref. [19] 
describe the expansion of the dependence in terms of 
Legendre polynomials and a convenient method of finding 
the m = component by integrating over the angle $. 



For A — 0, the e 1 term 



is 



1 = J < r ° ~ 3 ™) v-V2 (1 _ C0S 0)-V2, (B21) 
p V 2r2(r Q -2m) A V ; 1 > 



From Eqs. (C3) and (D7) in [19], this gives 



B4, = \ j 0m2 ° m — — ^ F 1/2 \/2. 



2r2(r G -2m) 



For A = 0, the e 1 term is 



(B22) 



' (r — 3m) 
2r2(r Q — 2m) 



1 (r Q + m) 

X ^oX 2 



X 1/2 (l-cos9) 



1/2 



(B23) 



From Eqs. (C3) and (D16) in [19] this gives 



-2V2 



3m) 



(2£-l)(2£ + 3) 4\/ 2r2(r Q - 2m) 



/2 



(r Q + m) 



3/2 



-2V2 



{2e-l){2£ + 3) 



(B24) 
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